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■ Abstract 
^ ■ 
O 
O _ 

The production of pions emitted backward in inelastic neutrino-nucleus 

o 

r~| . interactions is analyzed within the impulse approximation in the framework 

^ . of the dual parton model. We focus on the Q^-dependence of the multiplicity 

P5 . of negative pions, normalized to the total cross section of the reaction u + 



X 



A ^ ^-\- X. The inclusion of planar (one-Reggeon exchange) and cylindrical 



^ , (one-Pomeron exchange) graphs leads to a multiplicity that decreases as 

increases, in agreement with recent measurements carried out at CERN by the 

NOMAD collaboration. A realistic treatment of the high momentum tail of 

the nucleon distribution in a nucleus also allows for a satisfactory description 

of the semi-inclusive spectrum of backward pions. 
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Over the past several years, backward hadron production in inelastic lepton-nucleus scat- 
tering has been extensively investigated, both experimentally and theoretically 
with the aim of pinning down the dominant reaction mechanism and extracting information 
on the underlying dynamics. 

The analysis carried out on refs. focused on nuclear effects ine + A^e'+p + X 
reaction, while in ref. |^, hereafter referred to as I, we studied the spectrum of pions emit- 
ted in the process u + A — > ^ + tt + X . The results in I show that the main contribution 
to this spectrum comes from scattering off nucleons carrying high momentum. The cal- 
culations were performed within the impulse approximation (lA) scheme, using the Quark 
Gluon String Model (QGSM) P,p!0| to describe the interaction vertex. This model is similar 



to the Dual Parton Model (DPM) developed in refs. |TT|JT^ , in that both approaches are 
aimed at implementing the dual topological unitarization scheme. The relativistic invariant 
semi-inclusive spectrum of pions emitted strictly backward was obtained including only the 
cylindrical (i.e. one-Pomeron exchange) graph, while the contribution of the planar (i.e. 
one-Reggeon exchange) diagram was neglected. 

In this paper we continue the investigation of semi- inclusive pion production, started 
in I. We extend our analysis to pions emitted in the whole backward hemisphere, with 
respect to the direction of the incoming beam, and use QGSM to consistently include both 
the cylindrical and planar graphs in the calculation of the spectrum integrated over the 
transverse pion momentum. Our results are compared to the new data from the NOMAD 
collaboration [Q], which measured the spectrum of vr emitted in the backward hemisphere 
in the reaction u + C /x + vr + X and the Q^-dependence of the pion multiplicity, (ut^), 
normalized to the total cross section of the reaction u + C fi + X. 

Our study focuses on the kinematical region of large pion momentum {p.^ > 0.3 GeV/c), 
where the effects of final state interactions (FSI) leading to pion absorption associated with 
production of baryon resonances are expected to be small |l3| , p!^ and the lA can be safely 
used. 

The relativistic invariant semi-inclusive spectrum of pions produced in the process i + 



A ^ i' + IT + X, in which the incoming lepton is scattered with energy E' into the sohd 
angle dfl, is defined as 

PM^,,.x = E.^^^-^^ , (1) 

where E,, and p^^ are the total energy and three momentum of the produced pion, respec- 
tively. Within the framework of lA the quantity appearing on the right hand side of the 
above equation can be rewritten in convolution form, in terms of the semi-inclusive spec- 
tra of pions produced in lepton-proton and lepton- neutron processes pep^e'nx and /)^„^£/^x, 
according to (see, e.g., refs. P,P|,P^) 



PeA^e'7Txix,Q ■,z,p^t) = dy d kt fAiy.h) 

J x,z<y 



^ \y y y J 



"in^l'-n-X { — —,Pnt h 



(2) 



A \y y y _ 

where Z and N are the number of protons and neutrons, A = Z + N and = — g^, q being 
the four-momentum tranferred by the lepton. The relativistic invariant variables x and z 
are defined as 0,0] 

m 2{Pa-p,) ' ^ m (Pa ■ P.) ' ^ ^ 

where Pa and p^, are the four-momenta of the nucleus and the initial neutrino, respectively, 
and Ma and m are the nucleus and nucleon masses. The above variables are best suited to 
analyze the data of ref. 0, in which the backward direction is defined with respect to the 
beam direction. Note that they do not coincide with x and z employed in I, whose definition 
can be recovered by replacing the four-momentum of the incoming neutrino, p^, with the 
four momentum transfer q in eq.(|^). 

The nucleon distribution function fA^y, h) can be written P] 

fA{y, kt) = I dko dk, S{k) y5 (y - MAitPl^\ (4) 
J \ m [r'A ■ Pu) ) 

where Sik) is the relativistic invariant function describing the nuclear vertex with an outgo- 
ing virtual nucleon. The function resulting from the transverse momentum integration 



of fA{y,kt) can be obtained by approximating S{k) with the nonrelativistic spectral func- 
tion P(k, E) yielding the probability of finding a nucleon with momentum k and removal 
energy E = m — in the target nucleus |T^. However, due to the limited range of mo- 



mentum and removal energy covered by nonrelativistic calculations of P(k, E) (typically 
|k| < kmax ~ 0.7 — 0.8 GeV/c and (m — k^) < 0.6 GeV. (see, e.g., ref. IjT^)) this procedure 
can only be used in the region y < yo ^ 1.7 — 1.85. An alternative approach to obtain 
fA{y) at larger y, based on the calculation of the overlap of the relativistic invariant phase- 
space available to quarks belonging to strongly correlated nucleons, has been proposed in 
ref. [0. A similar procedure has been also used in ref. |[13| to obtain the quark distribution 
in deuteron at large y. 

The analysis of backward pion production requires the full nucleon distribution function 
given by eq.(^. We assume that it can be written in the factorized form 

fA{y,kt) = fA{y)9A{kt) , (5) 
with the function qa chosen of the Gaussian form 

where (fcj ) is the average value of the squared nucleon transverse momentum. 

The second ingredient entering the calculation of the spectrum of eq. @ is the elementary 
semi-inclusive spectrum. It can be calculated within the QGSM, based on the 1/A^ expan- 
sion, being the number of fiavors or colors, developed in |[T^. According to ref. 0, the 



elementary production process can be described in terms of planar and cylindrical graphs 
in the s-channel, as shown in fig. ^ for the case of neutrino interactions associated with the 
exchange of a IV-boson. The planar graph of fig. 1 (a) describes neutrino scattering off a 
valence quark. According to MM its contribution to the process u + N-^fi + Tr + X reads 



Fp {xi,Q'^;zi,pit) = zi(f)i{xi,Q^) 
with 



3 Vl-xi / 3 Vl-xi 



(7) 



Mxu Q') = Pfrr^d^ixu Q') , (8) 

71 I- xi + 

where G is the Fermi couphng constant, E is the energy of the incoming u, mw is the W- 
boson mass and xi = x/y = Q'^/2{k ■ p^) is the Bjorken variable associated with a nucleon 
bound in a nucleus. In eq.(|^), is the d-quark distribution inside a proton or a neutron, 
Duu^T^ and Dud-^-w are the fragmentation functions of uu or ud diquarks into positive or 
negative pions, pu = Pt — {z/y)kt is the transverse momentum of a pion produced on a 
nucleon carrying transverse momentum kt and Zi = z/y. 

According to refs. [p|,p!0| the planar graph of fig.l (a) corresponds to one-Reggeon ex- 
change in the t-channel, leading to the asymptotic behavior ly^^*^*^"* ^, where Wx is the 
squared invariant mass of the undetected debris and aij(O) = 1/2 is the Reggeon intercept. 
Hence, from 

Wx = {k + qf = k^ + Q'^^^Z^, (9) 

Xi 

where Xi = Q'^/2{k ■ q), it follows that the Regge behavior of the planar graph of fig.l (a) 
at moderate and large is given by 



The fact that the Q^-dependence of the graph of fig. 1 (a) is determined mainly by this 
Regge asymptotic is a consequence of the weak Q^-dependence exhibited by the calculated 



The second contribution to the spectrum, coming from the cylindrical graph of fig. 1 (b) 



can also be obtained within the approach of ref. |T0[. According to I, it can be written the 
form 

F^{xi,Q'^;zi,pu) = zi 02(Q^) [Li{xi,Q'^;zi,Pu) + Li{xi,Q'^; Zi,pu) 

+ L2{xi,Q'^;zi,pit) + L2ixi,Q^;zi,pit) , (11) 

with 



(12) 
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+ fdd{y, kt] Q^)Ddd-*7r f —,Pit - —kt 

\y y . 



Li = ds{xi,Q'^)Du^^{zi,pit) 



(14) 
(15) 



and 



1-^1 dy 
21 y 



d% dsiy,kt;Q')D 



-,Plt 



(16) 



y y J 

In the above equations Uy is the distribution of the valence w-quark, fud and fdd are the 
distributions of uu-, ud- and dd-diquarks inside the nucleon, Du-^-k, Dd-^n are the fragmenta- 
tion functions of the u- and (i-quark into a pion and Ddd^n is the fragmentation function of 
the (id-diquark dd into a pion. The exphcit expression of the quark and diquark distributions 
and fragmentation functions, obtained within the approach of ref. |10[, are given in I. 

According to refs. P,p!0| the diquark distribution fqqiy) coincides with the distribution of 
the corresponding valence quark evaluated at 1—y, i.e. with q^il—y)^ if qv{y) is normalized to 
unity. The main contribution to the pion spectrum is coming from target fragmentation, i.e. 
the kinematical region corresponding to large values of Zi = z/y. Within the approaches 
of refs. P,p!0| and |Tl],[12| the main contribution of sea quarks to the hadron spectrum is 
associated with n-Pomeron exchange processes (with n > 2), which are not taken into 
acount in this analysis as they provide a vanishingly small contribution at large y |T^-|T2[. 
We only include sea quarks produced from gluon decay, as shown in fig. 1 (b). In our 
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approach the gluon is seen as a nucleoli constituent, in addition to the three valence quarks. 
Therefore, we have used an average value of the gluon fraction < Ug > and normalized the 
valence quark distribution in a nucleon to 1 — (y^), the value {ug) ~ 0.15 — 0.2 being taken 
from ref. [0. This procedure amounts to setting fqqiy) = gi,(l — {i/g) — y). Our calculations 
show that the contribution of fragmentation of sea quarks and antiquarks (see the cylindrical 
graph of fig. 1 (b)), described by the terms Li and L2, is much smaller than that coming 
from fragmentation of valence quarks and diquarks, described by the terms Li and L2 of 
eq.dni). 

We assume a factorized form, similar to that of eq.(^, for the quark and diquark distri- 
butions and fragmentation functions. The dependence of the quark distribution is again 
chosen of the Gaussian form 

where (fc^j) is the average value of the squared quark transverse momentum. For the frag- 
mentation function we have used 



gqiqq)^H{kt) = I e-^'" . (18) 



The cylindrical graph of fig. 1 (b) corresponds to one-Pomeron exchange in the t-channel, 
whose asymptotic behavior is ly^^^^-*"^, where for the supercritical Pomeron P,pil| the value 



of the exponent is given by A = ap(0) — 1 ~ 0.1 — 0.15 |T^|T^. Hence, the Pomeron 
asymptotic of the cylindrical graph of fig. 1 (b) turns out to be {k"^ + (5^(1 — xi)/xi)^. The 
Q^-dependence of the cylindrical graph is dominated by the supercritical Pomeron behavior, 
as of eq. fpUD depends weakly upon . 

In conclusion, the relativistic invariant spectrum of pions produced in v + N /i + vr + X 
processes can be written in the form 

_ zd'^ J, , J {k^ + Q\l-x,)/x,) \-"' 
PuN^,nX = = Fp (xi, Q ; z,,p^t) 

(IXiCLZiCipit \ Sq / 



where sq = 1 (GeV/c)^ is a parameter usually introduced in Regge theory in order to get the 
correct dimensions. Substituting PuN^fj.TrX of eq.(|19D and the nucleon distribution fA{y,kt) 
into eq.(0) one can calculate the relativistic invariant spectrum of pions produced in the 
reaction u + A—^fi + ir + X. 

The multiplicity of pions normalized to the cross section of the process u + A ^ fj, + X , 
a, defined as 

{n-jr) = = - / dx — / d p^t PuA~^^l7^x , (20) 

can also be split into two parts, corresponding to the planar, or one-Reggeon exchange, 
diagram (fig. 1 (a)) and cylindrical, or one-Pomeron exchange diagram (fig. 1 (b)). While 
at fixed Xi the first term decreases as increases, the second increases with a slope dictated 
by the value of A. The Q^-dependence of (n^) resulting from our approach is shown in fig. 
H, together with experimental data taken from ref. 0. 

The NOMAD collaboration carried out a study of inelastic v-C interactions in which 
the negative pions emitted in the whole backward hemisphere, with respect to the incoming 
neutrino beam, were detected 0. The data of fig. show the multiplicity of negative pions 
carrying momenta in the range 0.35 < < 0.8 GeV/c, measured in a kinematical setup in 
which Wx increases as increases. Theoretical calculations have been performed applying 
the same kinematical conditions. 

The results obtained including both diagrams is represented by the solid line, while the 
dashed and dash-dot lines correspond to the separated contributions of the cylindrical graph 
of fig. 1 (b) and the planar graph of fig. 1 (a), respectively. 

The planar diagram provides the main contribution at small Q^, while the cylindrcal 
one dominates at > 10 (GeV/c)^. The decrease of the planar graph contribution to 
< n^r > and the increase of the cylindrical graph contribution with increasing are both 
consequences of the kinematical setup of ref. [Q], in which incresing leads to an increase of 
the squared invariant mass Wx (see eqs.(|^) and (|T^). It clearly appears that both diagrams 
have to be included to explain the observed Q^-dependence. The fact that the solid line 
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lies somewhat below the experimental data at > 2 (GeV/c)^ is likely to be ascribed to 
the contribution of secondary rescattering effects, which are not taken into account in our 
approach. 

The dotted line has been obtained setting P(k, £") = in the domain of large energy and 
large momentum, not covered by the nonrelativistic calculations of ref. ||TB[. Comparison 



between the solid and dotted line shows that the dominant contribution to (n^) comes from 
the high momentum tail of the nucleon distribution. 

Numerical calculations have been carried out using values of the average squared trans- 
verse momentum, entering eq.(§), in the range {kf) = 0.12 — 0.14 (GeV/c)^. These values 
correspond to an average nuleon momentum ~ 0.4 GeV/c, which seems to be reasonable for 
processes dominated by the high momentum tail of the nucleon distribution. The associated 
ambiguity in the results is always within the experimental errors. 

The other two parameters entering our calculations are the average quark and diquark 
transverse momentum and the slope of the Gaussian describing the fct-dependence of the 
fragmentation functions. Their values have been taken from ref. [^, where a modified 
QGSM, explicitly including the transverse motion of quarks and diquarks, has been devel- 
oped. 

Besides the dependence of the negative pion multiplicity, the NOMAD collaboration 
measured the semi-inclusive spectrum of negative pions, defined as 

27r^ da ^ 1 E^do_ 
a ""d^p^r o dpi ' 

in the kinematical region corresponding to pions emitted in the backward hemisphere with 
pI > 0.05 (GeV/c)^ [Q. In fig. |^ we compare the p^-dependence of the experimental 
spectrum to that resulting from our approach. It clearly appears that the contribution of 
the cylindrical graph dominates the spectrum and that the inclusion of the high momentum 
tail of the nucleon distribution, corresponding to p > 0.4 GeV/c, is needed to describe the 
data. 

In this paper we have analyzed inelastic neutrino-nucleus processes in which negative 



pions are emitted in the backward emishpere. The main conclusions of our work, concerning 
both the reaction mechanism and role of nuclear structure, can be summarized as follows. 
The dominant contribution to the reaction comes from target fragmentation, as shown by the 
calculation of the two topological graphs of fig. |1] within the dual topological unitarization 
scheme. Two kinds of the elementary processes have to be included: neutrino scattering 
off valence quarks, corresponding to one-Reggeon exchange in the t-channel (planar graph, 
fig. 1 (a)), and neutrino scattering off sea quarks, corresponding to one-Pomeron exchange 
in the t-channel (cylindrical graph, fig. 1 (b)). In the kinematical setup of ref. |^ the 
contribution of the planar graph to the analyzed pion multiplicity decreases as \JiJq^ as 
and Wx increase, while the contribution of the cylindrical graph increases as Q^^, with 
A in the range 0.1 — 0.15, as and Wx increase. The resulting Q^-dependence of the 
the normalized multiplicity of pions emitted in the backward hemisphere in the reaction 
V + C /i^ + TT" + X, turns out to be in fair agreement with the experimental data of 
ref. 0. As for the role of nuclear structure, the dominant contribution is coming from the 
high momentum tail of nucleon distribution, which can be described in terms of overlaps 
of distributions of three-quark colorless objects [Q. Including both graphs of fig. our 
approach also provides a satisfactory description of the p^-dependence of the measured pion 
spectrum. 

We are indebted to A.B.Kaidalov and M.Veltri for many helpful discussions. This work 
was partly supported by the U.S. Department of Energy. 
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FIGURES 





FIG. 1. Planar (a) and cylindrical (b) graphs contributing to the reaction u + A^ iJ, + h + X. 
Diagrams (a) and (b) describe processes in which the incoming neutrino interacts with a valence 
quark or a sea qq pair, respectively. 
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FIG. 2. Q^-dependence of the normalized multiplicity of pions produced in the backward 
hemisphere in the v + C — > /U^ + vr^ + X reaction. The solid line shows the results of the full 
calculation, including both graphs of fig. ||. The dashed and dash-dot lines correspond to the 
separated contributions of the cylindrical graph of fig. 1 (b) and the planar graph of fig. 1 (a), 
respectively. The dots show the results obtained setting P(k, E) = in the high energy-momentum 
domain, not covered by the nonrelativistic calculation of ref. The experimental data are taken 
from ref. IQ. 
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FIG. 3. p^-dependence of the spectrum of backward pions produced in the reaction 
1/ + C ^ //^ + TT^ + X. The solid curve corresponds to the fuh calculation whereas the dashed line 
has been obtained including the cylindrical graph of fig. 1 (b) only. The dash-dot line shows the 
results obtained setting ^(k, = in the high energy-momentum domain, not covered by the 
nonrelativistic calculation of ref. |jl6|. The experimental data are taken from ref. |^]. 
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